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We report results ol a study of the B 9 oscillation frequency using a large sample of B 9 semilep- 
tonic decays corresponding to approximately 1 fb _1 of integrated luminosity collected by the D0 
experiment at the Fermilab Tevatron Collider in 2002-2006. The amplitude method gives a lower 
limit on the B 9 oscillation frequency at 14.8 ps~ at the 95% C.L. At Am s = 19 ps _1 , the amplitude 
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deviates from the hypothesis A = (A — 1) by 2.5 (1.6) standard deviations, corresponding to a 
two-sided C.L. of 1% (10%). A likelihood scan over the oscillation frequency, Am,, gives a most 
probable value of 19 ps _1 and a range of 17 < Am s < 21 ps _1 at the 90% C.L., assuming Gaussian 
uncertainties. This is the first direct two-sided bound measured by a single experiment. If Am s lies 
above 22 ps _1 , then the probability that it would produce a likelihood minimum similar to the one 
observed in the interval 16 < Am B < 22 ps _1 is (5.0 ± 0.3)%. 

PACS numbers: 12.15.Ff, 12.15.Hh, 13.20.He, 14.40.Nd 



Measurements of flavor oscillations in the B® and B® 
systems provide important constraints on the CKM uni- 
tarity triangle and the source of CP violation in the stan- 
dard model (SM) [M. The phenomenon of B® oscillations 
is well established |2j, with a precisely measured oscilla- 
tion frequency Anid- In the SM, this parameter is pro- 
portional to the combination |V^Vtd| 2 of CKM matrix el- 
ements. Since the matrix element V ts is larger than V t d, 
the expected frequency Am s is higher. As a result, B® 
oscillations have not been observed by any previous ex- 
periment and the current 95% C.L. lower limit on Am s 
is 16.6 ps _1 0- A measurement of Am s would yield 
the ratio |Vt s /Vtd|i which has a smaller uncertainty than 
|Vtd| alone due to the cancellation of certain theory un- 
certainties. If the SM is correct, and if current limits 
on B® oscillations are not included, then global fits to 
the unitarity triangle favor Am s = 20.9^42 ps _1 Q or 
Am s = 21.2 ±3.2 ps" 1 0. 

In this Letter, we present a study of B®-Bg oscillations 
carried out using semileptonic B® — > /j, + D~X decays |5j 
collected by the D0 experiment at Fermilab in pp col- 
lisions at y/s = 1.96 TeV. In the B®-B® system there 
are two mass eigenstates, the heavier (lighter) one hav- 
ing mass Mh (Ml) and decay width Th (J^l)- Denoting 
Am s = M H - M L , AT S =T L - T H , T s = (T L + T H )/2, 
the time-dependent probability P that an initial B® de- 
cays at time t as B° -> n+X (P nos ) or B® -> yrX (P osc ) 
is given by p nos / osc = e~ Fst (l ± cos Am s t)/2, assuming 
that Ar s /T s is small and neglecting CP violation. Flavor 
tagging a b (b) on the opposite side to the signal meson 
establishes the signal meson as a B s ° a t time t = 0. 

The D0 detector is described in detail elsewhere 0. 
Charged particles are reconstructed using the central 
tracking system which consists of a silicon microstrip 
tracker (SMT) and a central fiber tracker (CFT), both 
located within a 2-T superconducting solenoidal mag- 
net. Electrons are identified by the preshower and liquid- 
argon/uranium calorimeter. Muons are identified by the 
muon system which consists of a layer of tracking de- 
tectors and scintillation trigger counters in front of 1.8- 
T iron toroids, followed by two similar layers after the 
toroids |jfl. 

No explicit trigger requirement was made, although 
most of the sample was collected with single muon trig- 
gers. The decay chain £?° — ► [i + D~X, D~ — > 07r~, 
4> — ► K + K~ was then reconstructed. The charged 
tracks were required to have signals in both the CFT 



and SMT. Muons were required to have transverse mo- 
mentum pt{h + ) > 2 GeV/c and momentum p(/i + ) > 
3 GeV/c, and to have measurements in at least two lay- 
ers of the muon system. All charged tracks in the event 
were clustered into jets ||, and the D~ candidate was re- 
constructed from three tracks found in the same jet as the 
reconstructed muon. Oppositely charged particles with 
Pt > 0.7 GeV/c were assigned the kaon mass and were 
required to have an invariant mass 1.004 < M(K + K~) < 
1.034 GeV/c 2 , consistent with that of a <p meson. The 
third track was required to have px > 0.5 GeV/c and 
charge opposite to that of the muon charge and was as- 
signed the pion mass. The three tracks were required to 
form a common Dj vertex using the algorithm described 
in Ref. 9]. To reduce combinatorial background, the Dj 
vertex was required to have a positive displacement in 
the transverse plane, relative to the pp collision point (or 
primary vertex, PV), with at least 4<r significance. The 
cosine of the angle between the Dj momentum and the 
direction from the PV to the D~ vertex was required to 
be greater than 0.9. The trajectories of the muon and 
D~ candidates were required to originate from a com- 
mon B® vertex, and the n + Dj system was required to 
have an invariant mass between 2.6 and 5.4 GeV/c 2 . 

To further improve B° s signal selection, a likelihood 
ratio method [lfj was utilized. Using M(K + K~tt) side- 
band (B) and sideband-subtracted signal (S) distribu- 
tions in the data, probability density functions (pdfs) 
were found for a number of discriminating variables: 
the helicity angle between the D~ and momenta in 
the (f> center-of-mass frame, the isolation of the n + D~ 
system, the x 2 °f the Dj vertex, the invariant masses 
M((j+Dj) and M(K+K~), and p T (K+K-). The final 
requirement on the combined selection likelihood ratio 
variable, y se \, was chosen to maximize the predicted ra- 
tio S/ VS 1 + B. The total number of D~ candidates after 
these requirements was iVtot = 26,710 ± 556(stat), as 
shown in Fig. QJa). 

The performance of the opposite-side flavor tagger 
(OST) [13 is characterized by the efficiency e = 
Nt&g/Ntot, where iVtag is the number of tagged B® 
mesons; tag purity T] s , defined as r) a — N COI /N tligl where 
N COI is the number of B® mesons with correct flavor 
identification; and the dilution T>, related to purity as 
T> = 2rj s — 1. Again, a likelihood ratio method was used. 
In the construction of the flavor discriminating variables 
x\,...,x n for each event, an object, either a lepton t 



5 



O6000 
(3 

o 

©4000 

w 

c 

>2000 
LU 



p» r\ n im it a 

- U\J Hun II / 

L 0) / 


\ 1 ID 







1200 



800 



400 



1.8 



1.9 
M 



(KK)it 



2.0 

[GeV] 



: D0 


Run II , 


\ 1 fb 1 


r w 












1.8 


1.9 


2.0 



FIG. 1: (K + K ~)tt~ invariant mass distribution for (a) the 
untagged B® sample, and (b) for candidates that have been 
flavor-tagged. The left and right peaks correspond to fj, + D~ 
and fi + D~ candidates, respectively. The curve is a result of 
fitting a signal plus background model to the data. 



(electron or muon) or a reconstructed secondary vertex 
(SV) , was denned to be on the opposite side from the B° s 
meson if it satisfied cosip(pi or sVjPb) < 0-8, where p B 
is the reconstructed three-momentum of the B® meson, 
and ip is the azimuthal angle about the beam axis. A 
lcpton jet charge was formed as Qj = YlitfPr/ J2iPri 
where all charged particles are summed, including the 
lepton, inside a cone of AR = (Aip) 2 + (Arj) 2 < 0.5 
centered on the lepton. The SV charge was denned as 
Qsv = E 1 (9Vl)°' 6 /E 1 (Pl)°' 6 > wh ere all charged parti- 
cles associated with the SV are summed, and p\ is the 
longitudinal momentum of track % with respect to the 
direction of the SV momentum. Finally, event charge 
is defined as Qev = J2i Q 1 Pt/ J2i Pt> where the sum is 
over all tracks with pt > 0.5 GeV/c outside a cone of 
AB, > 1.5 centered on the B° s direction. The pdf of each 
discriminating variable was found for b and b quarks us- 
ing a large data sample of B + — > n + vD° events where 
the initial state is known from the charge of the decay 
muon. 

For an initial b (b) quark, the pdf for a given variable 
Xi is denoted fi(xi) (/^(xj)), and the combined tagging 
variable is defined as d tag 

The variable dtag varies between 



(1 — z)/(l + z), where z 

nti(fK*i)/fi(.*i)) 

— 1 and 1. An event with d tag > (< 0) is tagged as a 6 
(b) quark. 

The OST purity was determined from large samples 
of B+ -> fi+D°X (non-oscillating) and B° d -> ^+ D*~ X 
(slowly oscillating) semileptonic candidates. An average 
value of eV 2 = [2.48 ± 0.21 (stat)+°;°g (syst)]% was ob- 
tained fll] |. The estimated event-by-event dilution as a 
function of |dtag| was determined by measuring T> in bins 
of Mtagl and parametrizing with a third-order polynomial 
for |(itag | < 0.6. For Idtagl > 0.6, V is fixed to 0.6. 

The OST was applied to the B® — ► /j, + D~X data sam- 
ple, yielding iV tag — 5601 ± 102 (stat) candidates having 
an identified initial state flavor, as shown in Fig. ^b). 
The tagging efficiency was (20.9 ± 0.7)%. 



After flavor tagging, the proper decay time of can- 
didates is needed; however, the undetected neutrino 
and other missing particles in the semileptonic B® de- 
cay prevent a precise determination of the meson's mo- 
mentum and Lorentz boost. This represents an impor- 
tant contribution to the smearing of the proper decay 
length in semileptonic decays, in addition to the res- 
olution effects. A correction factor K was estimated 
from a Monte Carlo (MC) simulation by finding the 
distribution of K = pt{h + D~)/pt(B) for a given de- 
cay channel in bins of M([i + D~). The proper decay 
length of each B® meson is then ct(B° s ) = ImK, where 
l M = M{B Q s )-{L T -p T {^ + D-))/{p T {^+D-)) 2 is the mea- 
sured visible proper decay length (VPDL) , Lt is the vec- 
tor from the PV to the B® decay vertex in the transverse 
plane and M{B° a ) = 5.3696 GeV/c 2 Q. 

All flavor-tagged events with 1.72 < M(K + K~-k~) < 
2.22 GeV/c 2 were used in an unbinncd fitting proce- 
dure. The likelihood, C, for an event to arise from a spe- 
cific source in the sample depends event-by-event on Im, 
its uncertainty 07 M , the invariant mass of the candidate 
M(K + K~ir~), the predicted dilution 2?(e? ta g), and the 
selection variable y sc \. The pdis for <Ji M , M(K + K~n~), 
T^(dtng) and y SG i were determined from data. Four sources 
were considered: the signal ii + D~{—t (j)ir~); the accom- 
panying peak due to /j, + D~(^ </>7r~); a small (less than 
1%) reflection due to n + D~(^ K + tt~it~), where the 
kaon mass is misassigned to one of the pions; and combi- 
natorial background. The total fractions of the first two 
categories were determined from the mass fit of Fig. ^b). 

The ^ + D~ signal sample is composed mostly of B® 
mesons with some contributions from B® and B + mesons. 
Contributions of b baryons to the sample were estimated 
to be small and were neglected. The data were divided 
into subsamples with and without oscillation as deter- 
mined by the OST. The distribution of the VPDL I for 
non-oscillated and oscillated cases was modeled appro- 
priately for each type of B meson, e.g., for £?°: 



„nos/osc 



p~~> - c (l,K, d tag ) 



(1) 



K 

CT B o 



exp(- 



Kl 

CT B o ' 



)[1 ± ©(dtag) cos(Am s • Kl/c)}/2. 



The world averages of t b o , t b + , and Arrid were used 
as inputs to the fit. The lifetime, t b o, was allowed to 
float in the fit. In the amplitude and likelihood scans 
described below, t b o was fixed to this fitted value, which 
agrees with expectations. 

The total VPDL pdf for the n + D~ signal is then the 
sum over all decay channels, including branching frac- 
tions, that yield the D~ mass peak. The B® — > n + D~ X 
signal modes (including D*q, and D s ^; and /i + 

originating from r+ decay) comprise (85.6 ± 3.3)% of 
our sample, as determined from a MC simulation which 
included the PYTHIA generator v6.2 [l2j interfaced 
with the EVTGEN decay package 0, followed by full 
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GEANT v3.15 modeling of the detector response 
and event reconstruction. Other backgrounds considered 
were decays via B° -► Df a) D~X and B°,B- -> DDj, 

followed by Df-. — ► /j + X, with a real Dj reconstructed in 
the peak and an associated real /i + . Another background 
taken into account occurs when the D~ meson originates 
from one b or c quark and the muon arises from another 
quark. This background peaks around the PV (peaking 
backgrounds). The uncertainty in each channel covers 
possible trigger efficiency biases. Translation from the 
true VPDL, I, to the measured Im for a given channel, is 
achieved by a convolution of the VPDL detector resolu- 
tion, of K factors over each normalized distribution, and 
by including the reconstruction efficiency as a function of 
VPDL. The lifetime-dependent efficiency was found for 
each channel using MC simulations and, as a cross check, 
the efficiency was also determined from the data by fixing 
Tgo and fitting for the functional form of the efficiency. 
The shape of the VPDL distribution for peaking back- 
grounds was found from MC simulation, and the fraction 
from this source was allowed to float in the fit. 

The VPDL uncertainty was determined from the ver- 
tex fit using track parameters and their uncertainties. To 
account for possible mismodeling of these uncertainties, 
resolution scale factors were introduced as determined by 
examining the pull distribution of the vertex positions of 
a sample of J/ip — > fJ. + ^ decays. Using these scale fac- 
tors, the convolving function for the VPDL resolution 
was the sum of two Gaussians with widths (fractions) of 
0.998cr; Af (72%) and 1.775o7 M (28%). A cross check was 
performed using a MC simulation with tracking errors 
tuned according to the procedure described in |15|. The 
7% variation of scale factors found in this cross check 
was used to estimate systematic uncertainties due to de- 
cay length resolution. 

Several contributions to the combinatorial back- 
grounds that have different VPDL distributions were con- 
sidered. True prompt background was modeled with a 
Gaussian function with a separate scale factor on the 
width; background due to fake vertices around the PV 
was modeled with another Gaussian function; and long- 
lived background was modeled with an exponential func- 
tion convoluted with the resolution, including a compo- 
nent oscillating with a frequency of Amj. The unbinned 
fit of the total tagged sample was used to determine the 
various fractions of signal and backgrounds and the back- 
ground VPDL parametrizations. 

Figure [2 shows the value of — Alog£ as a function of 
Am s , indicating a favored value of 19 ps -1 , while vari- 
ation of — log £ from the minimum indicates an oscil- 
lation frequency of 17 < Am s < 21 ps^ 1 at the 90% 
C.L. The uncertainties are approximately Gaussian in- 
side this interval. The plateau of the likelihood curve 
shows the region where we do not have sufficient resolu- 
tion to measure an oscillation, and if the true value of 
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FIG. 2: Value of — A log C as a function of Am s . Star symbols 
do not include systematic uncertainties, and the shaded band 
represents the envelope of all log C scan curves due to different 
systematic uncertainties. 

Am s > 22 ps -1 , our measured confidence interval does 
not make any statement about the frequency. Using 100 
parametrized MC samples with similar statistics, VPDL 
resolution, overall tagging performance, and sample com- 
position of the data sample, it was determined that for a 
true value of Am s = 19 ps -1 , the probability was 15% for 
measuring a value in the range 16 < Am s < 22 ps -1 with 
a — Alog£ lower by at least 1.9 than the corresponding 
value at Am s = 25 ps -1 . 

The amplitude method [l6J was also used. Equation Q] 
was modified to include the oscillation amplitude A as an 
additional coefficient on the cos(Ato s • Kl/c) term. The 
unbinned fit was repeated for fixed input values of Am s 
and the fitted value of A and its uncertainty o a found 
for each step, as shown in Fig. [3J At Am s — 19 ps -1 
the measured data point deviates from the hypothesis 
A = (A = 1) by 2.5 (1.6) standard deviations, cor- 
responding to a two-sided C.L. of 1% (10%), and is in 
agreement with the likelihood results. In the presence 
of a signal, however, it is more difficult to define a con- 
fidence interval using the amplitude than by examining 
the —A log C curve. Since, on average, these two meth- 
ods give the same results, we chose to quantify our Am s 
interval using the likelihood curve. 

Systematic uncertainties were addressed by varying in- 
puts, cut requirements, branching ratios, and pdf model- 
ing. The branching ratios were varied within known un- 
certainties Q] and large variations were taken for those 
not yet measured. The if -factor distributions were var- 
ied within uncertainties, using measured (or smoothed) 
instead of generated momenta in the MC simulation. The 
fractions of peaking and combinatorial backgrounds were 
varied within uncertainties. Uncertainties in the reflec- 
tion contribution were considered. The functional form 
to determine the dilution T>(d tag ) was varied. The life- 
time tbo was fixed to its world average value, and Ar s 
was allowed to be non-zero. The scale factors on the sig- 
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FIG. 3: -Bs oscillation amplitude as a function of oscillation 
frequency, Am a . The solid line shows the A = 1 axis for 
reference. The dashed line shows the expected limit including 
both statistical and systematic uncertainties. 



nal and background resolutions were varied within un- 
certainties, and typically generated the largest system- 
atic uncertainty in the region of interest. A separate 
scan of — Alog£ was taken for each variation, and the 
envelope of all such curves is indicated as the band in 
Fig. [21 The same systematic uncertainties were consid- 
ered for the amplitude method using the procedure of 
Rcf. 0, and, when added in quadrature with the sta- 
tistical uncertainties, represent a small effect, as shown 
in Fig. [3J Taking these systematic uncertainties into 
account, we obtain from the amplitude method an ex- 
pected limit of 14.1 ps -1 and an observed lower limit of 
Am s > 14.8 ps -1 at the 95% C.L., consistent with the 
likelihood scan. 

The probability that B®-B® oscillations with the true 
value of Am s > 22 ps" 1 would give a — Alog£ mini- 
mum in the range 16 < Am s < 22 ps -1 with a depth 
of more than 1.7 with respect to the — Alog£ value at 
Am s = 25 ps -1 , corresponding to our observation includ- 
ing systematic uncertainties, was found to be (5.0±0.3)%. 
This range of Am s was chosen to encompass the world 
average lower limit and the edge of our sensitive region. 
To determine this probability, an ensemble test using the 
data sample was performed by randomly assigning a fla- 
vor to each candidate while retaining all its other infor- 
mation, effectively simulating a B® oscillation with an 
infinite frequency. Similar probabilities were found using 
ensembles of parametrized MC events. 

In summary, a study of B®-B® oscillations was per- 
formed using B® — ► fi + D~X decays, where D~ — > 4>ir~ 
and cj> — > K + K~, an opposite-side flavor tagging algo- 
rithm, and an unbinned likelihood fit. The amplitude 
method gives an expected limit of 14.1 ps -1 and an ob- 
served lower limit of Am s > 14.8 ps -1 at the 95% C.L. 
At Am s = 19 ps — the amplitude method yields a re- 
sult that deviates from the hypothesis A = (A = 1) 



by 2.5 (1.6) standard deviations, corresponding to a two- 
sided C.L. of 1% (10%). The likelihood curve is well 
behaved near a preferred value of 19 ps -1 with a 90% 
C.L. interval of 17 < Am s < 21 ps -1 , assuming Gaus- 
sian uncertainties. The lower edge of the confidence level 
interval is near the world average 95% C.L. lower limit 
Am s > 16.6 ps~ x 0]. Ensemble tests indicate that if 
Am s lies above the sensitive region, i.e., above approxi- 
mately 22 ps -1 , there is a (5.0±0.3)% probability that it 
would produce a likelihood minimum similar to the one 
observed in the interval 16 < Am s < 22 ps -1 . This is 
the first report of a direct two-sided bound measured by 
a single experiment on the B® oscillation frequency. 

We thank the staffs at Fermilab and collaborating in- 
stitutions, and acknowledge support from the DOE and 
NSF (USA); CEA and CNRS/IN2P3 (France); FASI, 
Rosatom and RFBR (Russia); CAPES, CNPq, FAPERJ, 
FAPESP and FUNDUNESP (Brazil); DAE and DST 
(India); Colciencias (Colombia); CONACyT (Mexico); 
KRF and KOSEF (Korea); CONICET and UBACyT 
(Argentina); FOM (The Netherlands); PPARC (United 
Kingdom); MSMT (Czech Republic); CRC Program, 
CFI, NSERC and WestGrid Project (Canada); BMBF 
and DFG (Germany); SFI (Ireland); The Swedish Re- 
search Council (Sweden); Research Corporation; Alexan- 
der von Humboldt Foundation; and the Marie Curie Pro- 
gram. 



[*] On leave from IEP SAS Kosice, Slovakia, 
[f] Visitor from Helsinki Institute of Physics, Helsinki, Fin- 
land. 

[1] S. Eidelman et al, Phys. Lett. B 592, 1 (2004). 

[2] E. Barberio et al. (Heavy Flavor Averaging Group), 

hep-ex/0603003 Note that we take h = c = 1, hence 

the units on Am,. 
[3] J. Charles et al. (CKMfitter Group), Eur. Phys. J. C41, 

1 (2005). 

[4] M. Bona et al. (UTfit Collaboration), J. High Energy 
Phys. 07(2005) 028. 

[5] Charge conjugate states are assumed throughout. 

[6] V. Abazov et al. (D0 Collaboration), physics/0507191 
[Nucl. Instrum. Methods Phys. Res. Sect. A (to be pub- 
lished)]. 

[7] V.M. Abazov et al, Nucl. Instrum. Methods Phys. Res. 

Sect. A 552, 372 (2005). 
[8] S. Catani et al., Phys. Lett. B 269, 432 (1991), 

"Durham" jets with the pr cut-off parameter set at 

15 GeV/c. 

[9] J. Abdallah et al. (DELPHI Collaboration), Eur. Phys. 

J. C 32, 185 (2004). 
[10] G. Borisov, Nucl. Instrum. Methods Phys. Res. Sect. A 

417, 384 (1998). 
[11] V. Abazov et al. (D0 Collaboration), Phys. Rev. D 

(to be published); D0 Note 5029, available from 

http: / /www- d0.fnal.gov/Run2P hysics/WWW| 

/results/prelim/B /B32/ . 



8 



[12] T. Sjostrand et al, Comput. Phys. Commun. 135, 238 
(2001). 

[13] D.J. Lange, Nucl. Instrum. Methods Phys. Res. Sect. A 

462, 152 (2001). 
[14] R. Brun and F. Carminati, CERN Program Library Long 

Writeup W5013 (unpublished). 



[15] G. Borisov and C. Mariotti, Nucl. Instrum. Methods 

Phys. Res. Sect. A 372, 181 (1996). 
[16] H.G. Moser and A. Roussarie, Nucl. Instrum. Methods 

Phys. Res. Sect. A 384, 491 (1997). 



